State Transfer in Highly Connected Networks and a Quantum Babinet Principle 
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The transfer of a quantum state between distant nodes in two-dimensional networks, is considered. 
The fidelity of state transfer is calculated as a function of the number of interactions in networks 
that are described by regular graphs. It is shown that perfect state transfer is achieved in a network 
of size N, whose structure is that of a -y-cross polytope graph, if N is a multiple of 4. The result 
is reminiscent of the Babinet principle of classical optics. A quantum Babinet principle is derived, 
which allows for the identification of complementary graphs leading to the same fidelity of state 
transfer, in analogy with complementary screens providing identical diffraction patterns. 



I. INTRODUCTION 

The dynamics of quantum many-body systems offers 
a rich variety of features. This quantum dynamics is 
often investigated in one-dimensional chains, which are 
amenable to exact analytical treatment in some cases 
1] and, in other cases, efficient numerical methods have 
been developed for their simulation [2J . For more general 
lattice structures, henceforth referred to as graphs, few 
analytical treatments are known. An important prob- 
lem that arises in this context is the interplay between 
the dynamics of quantum many-body systems and the 
properties of the underlying graph, which determines the 
interaction structure of the many-body system. A vari- 
ety of interesting phenomena, two examples of which are 
perfect state transfer [J Q and thepossibility of deciding 
the graph isomorphism problem [5j, have recently been 
explored in such complex quantum systems. 

In the case of graphs with uniform nearest neighbour 
coupling, perfect state transfer (PST) has so far been 
proven possible only with rings of N = 4 spins, chains 
of N = 2 or N = 3 spins and with Cartesian products 
of such graphs, the so-called one-link and two-link hy- 
percubes [J d, 0- For larger networks, it appears that 
increasing the number of spins and the degree of the un- 
derlying graph tends to compromise the transmission of 
quantum information [H, 0, In the static case it 

has been shown that higher connectivity and associated 
monogamy constraints frustrate the system and affect its 
quantum correlations [Til ]. On the other hand, if natural 
interactions are abandoned in favor of particular coupling 
schemes, in which only nearest neighbors interact and the 
interaction strength depends on their position relative to 
a fixed point, then perfect state transfer is possible in 
spin chains with large N 0, [H, 0, • 

In this paper we investigate whether it is possible to 
transfer perfectly a quantum state between two distant 
nodes of a two-dimensional spin network, in which the 
interactions between spins are both permanent and ho- 
mogeneous. We show that PST can be achieved in such 
a two-dimensional highly-connected network of arbitrary 



size, N. This is possible with a unique regular configu- 
ration, namely a two-dimensional graph of the -^-cross 
polytope [HI], which is dual to the hypercube in ^ di- 
mensions and isomorphic to a type of circulant graph 
[l6|. It turns out that these findings lead to a natu- 
ral quantum generalization of a well-known principle in 
classical optics. Therefore the plan of the paper is the 
following: in Sec. [TXJ we introduce a general spin model, 
whose defining characteristic is that it preserves the total 
number of excitations in the network; then in Sec. IIIII wc 
present numerical calculations, which reveal the special 
properties of Y _cross polytope graphs; and in Sec. HVl 
we provide analytical results that support our numeri- 
cal findings and prove the main result of the paper. In 
Sec. El based on the quantum state transfer properties of 
complementary graphs, we derive a quantum version of 
the Babinet principle from classical optics. Basic results 
of Monte-Carlo simulations on the influence of static dis- 
order on the system are presented in Sec. IVIl Concluding 
in Sec. IVII1 we discuss our results. 



II. EXCITATION-PRESERVING QUANTUM 
NETWORK 



We begin by considering N spins- \ situated along a 
circle, as shown in Fig. [TJ It is understood that if two 
spins are interacting, a line is drawn between them. The 
result is a graph Q — (V,E); the vertices V{Q) represent 
the spin sites and the edges E(Q) represent pairwise in- 
teractions. The necessary information about the graph 
Q is contained in its adjacency matrix, A(Q), whose ele- 
ments are given by A^j = 1 if € E(Q) and are zero 
otherwise. We consider Hamiltonians of the form (H = 1) 



JY 



H 



1, (1) 



fc=l 



where (cr^T) are the raising and lowering operators 
for site k, ujk is the local site excitation energy and J k ,i 
denotes the hopping rate of an excitation between the 
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FIG. 1: (Color online) A circulant spin network can be used to 
transfer quantum states from A to B. In a network of N — 6 
spins there are three possible configurations with connectivity 
(a) C = 1, (b) C = 2, and (c) C = 3, as shown. Network (b) 
is a 3-cross polytope graph (CPG). 



sites k and I. The dynamics in this system preserves 
the total excitation number, defined by N'=^2k=i a k a k • 
During dynamical evolution the state of the network, 
|*(t)) = exp(-iHt)\^ ), where |* ) is the initial state, 
always remains in the same excitation sector because 
[H,Af] — 0. In what follows we restrict our attention 
to the single-excitation sector, for simplicity. In this sub- 
space the Hamiltonian of the system is equal to the ad- 
jacency matrix of the underlying graph, Ti = A(Q), pro- 
vided that the spin-spin interactions are homogeneous. 
Deviations due to engineering errors in the interactions 
are also examined later on. The network is prepared in 
the state 



\^o(j)) = \j) := •>,(>,••• I..-" 



(2) 



where only spin j is excited. The propagation of an arbi- 
trary state a\0j) + where \a\ 2 + \f3\ 2 — 1, is equiv- 
alent to the propagation of the state (since the +1 
eigenstate of Zj, \0j), does not evolve under TL). The aim 
is to transfer the excitation from j to N/2 + j, that is, 
to the vertex that is diametrically opposite from j across 
the ring - hence we initially consider that N is even. The 
state transfer is quantified by the fidelity 



F(t) := K* (N/2+j)|exp(-iWt)|*oa)>l 



(3) 



Perfect state transfer is achieved at a certain time to if 
and only if F(t a ) = 1. 



III. QUANTUM STATE TRANSFER AND 
CONNECTIVITY 



We now ask, 'How is the fidelity of quantum state 
transfer influenced by the connectivity of a network?' 
The connectivity, C(Q), is defined here as the number of 
edges that are incident on a vertex, counting only within 
the half-disc defined by that vertex and the opposite one 
(i.e., it is half the degree of the graph). The extreme cases 
are those of a ring (C = 1) and a fully-connected network 
(C = N/2), but in general we have C = 1, 2, • • • , N/2 
(see Fig. [1]). Before we analyze this question analytically 
we calculate numerically the fidelity F(t) of Eq. © for 



t £ [0, At], given the number of spins N and the connec- 
tivity C. The maximum fidelity, max(i 7 At)' is then de- 
termined for the interval At. It is assumed that tut = 
and Jk,i = 1 in the Hamiltonian of Eq. (fl]). Therefore 
the Hamiltonian of the network is equal to the adjacency 
matrix of the underlying graph structure. Under these 
conditions it is observed in Fig. [2a) that the fidelity is 
a non-monotonic and rather complicated function of the 
connectivity. However, it displays remarkable behavior 
for C = ^ — 1, which corresponds to a 2£;-cross polytope 
graph (CPG) with N = 4fc spins, where k is a positive 
integer. In this case, PST is achieved at to = § + n7r , 
i.e., F(^ + mr) = 1, where n > is an integer. For k = 1 
we recover the known result [l7| for a ring with N = 4. 



In Ref. [l6( it was shown that circulant graphs of odd 
order do not allow perfect state transfer (so our choice of 
even N is justified) and, moreover, it was left as an open 
question if there exist circulant graphs of even order with 
N > 4 that support PST. Our results show that such 
graphs do indeed exist: the 2fc-CPG is isomorphic to the 
circulant graph Ci4fc(l, 2, . . . , 2k — 1). In these networks 
every spin interacts with every other spin, except for one 
(e.g., see Fig. [IJb) for an example). The appropriate 
choice of At is made by comparing trial values with the 
occurrence time of the first peak in the evolution of the 
fidelity for a spin ring (this evolution is shown in Fig. 
[5fb)). In Fig. [21(c) we show the evolution of the fidelity 
for a network with connectivity C — y — 1 = 99. It is 
seen that the fidelity becomes equal to 1 at to = k/2. 



(a) N=200, J=1,At=100 




FIG. 2: (Color online) (a) Maximum fidelity in the interval 
[0, At = 100] against connectivity for a network of size JV = 
200 and homogeneous interactions, (b) Fidelity against time 
for the simple ring network, (c) Fidelity against time for the 
100-cross polytope graph network. 
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IV. ANALYTICAL RESULTS 

In this section we analyze the perfect (for N = 4k) or 
near perfect state transfer (for TV = 4fc + 2) in the config- 
urations introduced previously. The Hamiltonian Hcpg 
of a cross polytope graph is that of Eq. (pj with ujf. = 
and J k) i = (1 - S l N +k )(l - S k j) = J lM . The Hamil- 
tonian of a fully-connected network, denoted as Hf Cl is 
that of Eq. (JTJ) with J k j = 1 — S k j = Ji k . The Hamil- 
tonian where only opposite pairs are connected, denoted 
as H pair , has J k j = S t N_ +k = J tjk for all k = 1, . 

We start by noting that we have 

HcPG = Hf c — Hp a i r , [Hfc, H pair ] = 0. 

Therefore, the evolution operator is 



AT 

2 



e -iH C PGt — e iH pair t e -iHf c t 



Due to the fact that 

N/2 



fc=l 



(4) 
(5) 

(6) 



only couples opposite sites, its dynamics is very simple. 
It is straightforward to obtain 



e iH pazr t = icost . 



N/2 ( N 

l E l fc )(y + fc l 

k=l ^ 



k)(k\ sint. 



In particular, after a time t = % + mr a starting state \k) 
will have been transformed to (— + k). Finally, 
to determine the dynamics of He pg we need to consider 
Hf c . The latter can be expressed as 



H 



fc 



N\+)(- 



1. 



where 



N 



(7) 



(8) 



We have 



e- lH ^ f = [\+)(+\e~ lNt + (1 - l+X+D] e l *. 

Therefore, a state \k) is mapped onto itself, up to a global 
phase, under e ~ lH f at when Nt = 2wk with k £ N. 

As a consequence, the dynamics due to Hcpg allows 
for PST if both Nt = 2-nk and t = | are satisfied for the 
same t. This implies the condition 



N = 4fc 



(9) 



and explains the possibility of PST in 2/c-cross polytopes. 
For N — Ak + 2 the analysis above immediately applies 
and shows that we do not have PST at t = ^. 



More generally, we can find the transfer fidelity for 



t = ^. Starting with \k) and using 
we find at t = \ the state 



„iH„ 



rt e -lH fc t\ 



-_L|+) e -^/2 + 

N 



N 
2 



e iH pair tl 



1 xfN 



Then the fidelity, |(f + k\e iH ^ irt e- iH f^\k) | 2 , is 



N \ N 



1 — cos • 



Ntt 



(10) 



For N — 4k we recover F = 1, while for TV = 4fc + 2 we 
find that F = (1 - 2/7V) 2 . Therefore, as TV — > oo, the 
fidelity approaches 1 and we obtain almost PST. 



V. QUANTUM BABINET PRINCIPLE 

These results provide a clear insight into the basic 
mechanisms that facilitate PST in these systems. The 
key realization is that a fully connected network in which 
some couplings J k j are removed, can behave similarly to 
an initially unconnected network which is supplemented 
with the very same J k j links. This result is in fact 
reminiscent of the Babinet principle of classical optics 
[HI], which is illustrated in Fig. [3J In our context of 
state transfer through connected networks, the situation 
is similar in the sense that 



e - lHc PGt e - 



-iHfat 



because Hcpg and H pair commute and also e ~ lHfct 
equals the identity at specific times t (in the opti- 
cal setting this is the situation when all incident light 
emerges unaffected). Of course in the quantum setting 



we have the added problem that e lH cpGt 



iH„ 



e -i(H C PG+H pair )t in genera l_ 

The analog of the Babinet theorem does hold however 
for much more general settings than just that of commut- 
ing Hcpg and H pa i r . Indeed, as before, let us assume 
that 



H f = NP, P=|+)(+| 
For a sequence of Hm that satisfies 

\\FH N (t-F) + {l-F)H N 



= 0(-==) 
WN 



we compare the dynamics of H^ and Hf — H 



IN 



(11) 



(12) 



H, 



in the limit of large N. The following argument is not 
fully rigorous in that it does not provide detailed error 
estimates but these may be provided in a more detailed 
analysis. 

Let us consider the dynamics under H c in an interac- 
tion picture with respect to Hf when this Hamiltonian 
becomes time-dependent, 

H cJ (t) = e iH i\H c - H f )e- iH f l = -e^^Hm^^ . 
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FIG. 3: (Color online) Illustration of Babinet's principle in 
an optical set up with Fraunhofer conditions, (a) An unob- 
structed plane wave is focused by a lens L and produces a 
diffraction pattern of amplitude A(r) on the screen S. (b) 
Diffraction patterns resulting from complementary screens s 
and s, whose opaque and transparent areas are swapped. At 
any point downstream from s and s, the sum of the two 
diffracted amplitudes, A s (f) + As(r), equals the amplitude 
diffracted from the unobstructed plane wave A(r). Away 
from the central spike, this amplitude is zero and therefore 
A a (f) = — As(f), which leads to Babinet's prediction of identi- 
cal diffracted light fields for complementary apertures. Com- 
plementary apertures play the role of complementary graphs 
describing quantum spin networks. Increasing r corresponds 
to increasing the number of nodes N. 

The corresponding time-evolution operator from t\ until 
<2 in the interaction picture will be denoted by t/j(ii — > 
£2)- Now we note that H c j(t) may contain rapidly os- 
cillating terms (those coupling subspace defined by P to 
the subspace defined by 1 — P) thanks to the action of 
e f*. These rapidly oscillating terms may be neglected 
for large N leading to a correction of order 1/N in the 
dynamics. Hence, we find 

H c ,i(t) -(1-F)H N (1-F) -FH N F (13) 
= -H + FH N (t -F) + {1 -F)H N F. 

As we had assumed earlier that \\FH N (t - P) + (1 - 
F)H N F\\ is of order 1/VN we find that H c j(t) is well 
approximated by — up to corrections that decrease 
with increasing N. Hence (7/(0 — > t) w e tHNt , and we 
find that 

e -iH c t = Q-iB/tjj^Q _> t ) s e -iH f t e iH N t_ ( W ) 

Now we consider the transfer fidelity from state 
\k) to l^ + fe)) a s an example. The amplitude 
(f + k\e- iH ^\k), using e"*^* = e -iH t t e iH N t^ ig found 
to be equal to 

(f + k\e iHNt \k) + (y + k\+)(+\e iHNt \k){e lNt - 1) 

<*(^ + k\e iH »*\k) 

where the difference decreases with increasing N. There- 
fore, the transition amplitudes according to the dy- 
namics under —Hn and H c are asymptotically (in N) 



equal. Note that for a real Hamiltonian Hpj we have 
((N + k\e lHNt \k))* = (f + k\e~ iHNt \k) so that 

|<^+fc|e-«^|fc>| - |(^ + fe| e -^*|fc)| (15) 

again with an error that decreases with increasing N. 
This is the quantum Babinet principle. 



VI. INFLUENCE OF DISORDER 

We provide here a brief analysis of realistic engineering 
errors in the interactions of a 2fc-CPG network in order 
to assess the robustness of a possible experimental im- 
plementation. We take into account two types of errors: 

(i) disorder in the interactions, and (ii) random breaking 
of interactions. For case (i) we assume that if p and q 
are interacting then the interaction strength can take any 
value in the interval [1 — 8, 1 + 8] , with equal probability. 
The amount of disorder is thus quantified by 8 € [0,1]. In 
case (ii) some interactions are randomly broken, that is, 
J pq vanishes for a fixed number of pairs (p, q) . The num- 
ber of broken interactions is B E [0, 1), given as a ratio 
to the total number of interactions in the network. The 
main results of Monte-Carlo simulations on ^-CPG net- 
works with N = 40, 80, 120, 200, 400 spins, are as follows. 
For type-(i) errors we find that disorder up to 8 = 0.02 
allows for almost PST in smaller networks (N < 100). In 
particular, the maximum fidelity F is greater than 0.99, 
on average, with a worst-case value of 0.98 in the case 
of N = 40; while for N > 100 the average maximum F 
is over 0.95 for disorder that is less than 2%. For type- 

(ii) errors we find that the random breaking of very few 
bonds, so that B < 0.001, still allows for very high qual- 
ity state transfer, where the maximum F is larger than 
0.95, on average. However, the value of the worst-case 
fidelity peak fluctuates considerably on individual cases, 
depending on the positions of the broken bonds. 

In this connection, the usefulness of the quantum Babi- 
net principle can be illustrated in the case of transport 
of excitations through noisy networks, a setting that has 
recently been introduced independently in [lj| and [2fj| . 
Initially all population resides in a given site and we eval- 
uate how much population may be transferred asymptot- 
ically to a selected target site. To this end, we let the 
target site be attached to a sink to which the population 
is transferred irreversibly. We want to analyze whether 
the presence of local dephasing can assist the excitation 
transfer. If the sink is attached to site ^ + 1, then the 
Babinet principle implies that the evolution is that of a 
system where only the opposite two sites are coupled, 
and we recover a situation for which it was proven in J2CM] 
that no dephasing enhanced transport is possible [211 ]. 
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VII. SUMMARY AND DISCUSSION 

We have shown that PST is achieved in a network of 
size N, whose structure is that of a Y~ cross polytope 
graph, if N is a multiple of 4. If N is even, but not a mul- 
tiple of 4, then almost PST is achieved for larger networks 
of this kind, so that F approaches 1 for N — ► oo. These 
results can be interpreted in terms of a quantum Babinet 
principle, which establishes the conditions required for 
having complementary graphs leading to same fidelity of 
state transfer, in analogy with the classical situation of 
obtaining identical diffraction patterns from complemen- 
tary screens. As shown in various examples, invoking 
Babinet 's principle alone can simplify the analysis of the 
performance of connected networks and therefore become 
a useful tool in tackling a variety of problems in quantum 
information theory. 
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